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Abstract 

The improved method of intermittent data analysis is proposed. It ex- 
ploits, in addition to the standard density moments, the information on the 
bin-bin correlations, observed in the data and expressed in terms of the den- 
sity correlators. The improving recovering corrections are implemented into 
the data analysis in the form of the recursive algorithm, and tested in the 
framework of multiplicative cascading models. 

The first signals on possible intermittent behaviour in high-energy multiparticle 
production M were found in the data of the single event recorded by the JACEE 
collaboration . The presence of large dynamical fluctuations manifesting a scaling 
behaviour was registered also afterwards in other accelerator experiments |J . Many 
different models [|J have been proposed since to explain the power-law rise of the 
multiparticle moments, described by the coefficients called the intermittency expo- 
nents. Some of the models suggested that the observed scaling may be the result of 
final state multiparticle cascading and the intermittent data represent the last 
stage of the cascade. In this approach the main problem lies in the extraction of the 
information on the previous cascade stages which are in some way encoded in the 
last stage data. The standard method of recovering the history of the cascade was 
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proposed and applied originally to the JACEE event data. Since that time it has 
become a standard tool of multiparticle data analysis [f|], especially in the event-by- 
event analysis |5j, || . However, the theoretical problem what is the interplay between 
the cascade recovered in the standard analysis and the true cascade which generated 
the data has been studied only recently 0, ||. It ended up with the proposal of 
introducing the new recovering corrections to the data analysis. 

In this paper we would like to summarize the results on improvement of the 
standard data analysis achieved by including the recovering corrections. Most of 
these results were derived in Ref. ||. In what follows we will concentrate on the 
technique of implementing the recovering corrections into the data analysis in the 
way which may be useful for experiment. 

The recovering corrections aim to improve the method of recovering the his- 
tory of the particle cascade, also called the rebinning M. The standard method 
reconstructs the cascade from the last step data, represented by a sample of M 
numbers : xf 1 ^ (i = 1,...,M). They describe e. g. the distribution of particle 
density into bins. For simplicity assume M = 2™, where n denotes the number of 
cascade steps. The recovered particle density yf l ® in the zth bin at the (n — k)th 
cascade step (k — 0, . . . , n — 1) takes then the form : 

i 2 k -l 

(n-k) _ J_ (n) /. s 

y « ofe • h 2 k xi+j- \ L ) 

Z 3=0 

The intermittency exponents are extracted from unnormalized reconstructed den- 
sity moments z q k ' rec [] : 

Z 3=0 

assuming that z^\ ec manifests a power law behaviour : 

(*) _ 2 k <, (3) 

q;rec. ) \ u / 

and the normalized intermittency exponent (p q - n orm. = <f>' q , where <p q;n orm. '■= — 
q<pi. The exponent 4>' q is estimated as a slope of the linear x 2 — fit applied to the 
points (k, log z q k l ec ) (logs = log 2 x). 

However, it was found and proved in 0, || that there exists a difference between 
the true density moments : 

z i=0 

1 Actually, the factorial moments are normally used (to reduce the statistical noise) but this 
does not concern us here 
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and the reconstructed ones (|3|) which influences the estimation of intermittency 
exponents. This difference may be expressed in the form of correcting factor p q (k) : 

which was called the recovering correction ||. 

The recovering corrections contain information on the specific process which 
generated the true cascade. In what follows we restrict ourselves to the corrections 
considered in the framework of multiplicative random cascading models || |T7j, |TTJ 
which are nowadays widely recognized in multiparticle data analysis ||. Similarly as 
in ||, we consider the multiplicative models with possible neighbour-node memory 
which generate a uniform distribution of particle density into bins. In the multi- 
plicative cascade the root of the cascade is set equal 1 : Xq ^ = 1. One generates the 
next stages of the cascade recursively, following the scheme : 

:= W 2 .x?\ (6) 

where W\ and W 2 are random variables of m model parameters aj, j = 1, . . . , m : 

W\ = ctj with probability p a . , 
W 2 = (ij with probability p a . , 

(7) 

with normalized probability weights Y^JLiPaj — 1- The distribution of particle 
density will be uniform if the following condition is fulfilled : 

p{W x = at, W 2 = aj) = p{W x = a h W 2 = a*), (8) 

where p(W\ = ai,W 2 = aj) denotes probability of choosing in W\ = a>i and 
W 2 = aj = l,...,m). Then unnormalized density moments fulfill the 
scaling relation : 
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(fc) ^ 2 fc '^, (9) 



and intermittency exponents <p q read : 

(j) q = \og(a q 1 p ai + ... + a q m p am ). (10) 



The popular models : a—, p— models 0,|TT[ and the (p + a)— model introduced in 
II are special cases of multiplicative rule @. 

It was proved in || that for any multiplicative process which obeys rules (^j), @ 
recovering corrections p q (k) fulfill the recurrence equation f\ : 

- \ q E ( ) ) V 3 ik - l)p q _j(k - l)(W{Wr 3 ) (11) 



2 A similar recurrence relation has been obtained in a different context in |li 
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with the initial conditions 



p g (0) = 1, 

Po(k) = 1. (12) 



We introduce a notation : 



(W(W l 2 ) = k hl . (13) 

Formula (|i~l~D implies that coefficients kjj are the only parameters of the multiplica- 
tive model needed for calculating the value of p q {k). Furthermore, it is not difficult 
to establish the values of kjj from the model. For either j = or I = they equal : 

/■-,.» luij ■>''■. (14) 

where 0/s are ordinary intermittency exponents @. To find the value of kjj for 
both j, / ^ we use the unnormalized density correlators 0, £|, [13 



c j,i ~ ofc-i X! ( x 2i' > ) • (15) 

The correlators and the density moments fulfill the relation (see ||) : 

C$1 = # + T 13 • ^, (16) 

which can be also rewritten as : 

logcg ) = (A;-l)0 i+i + log% i . (17) 

Both relations (16), imply that we may derive fc^y in a straigthforward way by 
calculating correlators and density moments from data, and applying to them the 
standard x 2— fit- 
Applying the standard method to the correlators at the previous cascade stages, 
we would expect to find the similar difference between the reconstructed correlators 
and the true ones, as it was observed for the density moments. It was proved in 
H that this difference may be expressed in terms of the same recovering correction 
p q (k) (|TT|) as for the density moments : 

c^i'rl. = c^i~ k) p j+ i(k). (18) 

Now we have all tools needed for implementation of recovering corrections into 
the multiplicative data analysis. To illustrate the problem we describe the improved 
estimation of intermittency exponents of the second rank. The recovering correction 
of the second rank derived from (ITTD reads : 



p 2 (k) = \{p 2 (k - l)fc 2 , + P2 (k - l)fc 0)2 + 2p\{k - l)fc M ) (19) 



with the initial condition p2{0) = 1. It was proved in H that : 



pi (k) = 2 k * 1 = (zi n) ) n (20) 

The parameters needed to calculate P2(k) are following : z[ n \ 2 (cf-(0)) an d 
Derivation of z\ is straigthforward, and to estimate the values of </> 2 and ki t i we 
propose the following recursive procedure. The primary values of (f>2 and k^i may 
be obtained in the standard way from CPD ? (P~^D • We substitute them to formula 
( pT)|) to derive the approximate form of correction P2{k). Now, applying again the 
approximate form of P2(k) to equations @, (|16D , one derives adjusted parameters 
4>2 and fei^ and compares them with the primary values. If the relative difference is 
large, one repeats the recursive adjusting till the parameters do not change within 
a given accuracy. 

One may generalize the above scheme for the intermittency exponents of any 
rank. Following ||, below we present the implementation algorithm which recur- 
sively adjusts the primary parameters <p q , kj 7 i (j + I = q, jl > 0) obtained after 
applying the standard method to the data : 

(INPUT) parameters <p ± , . . . , <f> q _i, kjj (j + I = 1, . . . , q — 1) obtained after 
applying the implementation algorithm for q — 1, 2, . . . , q — 1 step- by-step 

(1) derive <f)' , k'j q _j (j = 1, . . . , q — 1) from data, using the standard method 
i. e. reconstruct the cascade using (|l|) and derive the parameters from relations : 

log 41c. = + (21) 



k 



(k) _ (fc-1) 

L"i /■ von «i 



rec. j+l', Tec. j,l 



■ K, (22) 



where k = l,...,n (cf. (|), ©), 

(2) derive <f> q;corr . , kj A -j. iCarr . (j = 1, - . - , g — 1) in the following substeps : 
(2.0) calculate p q (k) from relation (cf. (PD ) : 



1 A, q 



2? ^ \ 7 



using 0' derived in step (1), and estimate <fi q . corr . from : 

log ~ hg(p q (k)) = (n-k)- <\> q , corr . + b, (24) 
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(2.1 ) calculate p q (k) from using (fi q - CO rr. (other parameters as after step (1)), 
and estimate corn from relation (cf. (|17|),([2"3"1)) : 



\°% C tl-,rec. ~ l °S{Pj+l{k)) = {n-k- + log kjj-corr., (25) 



(2.q-l) calculate p q {k) from (|23|), using all previously derived parameters cf) q . corr , 
kj, q ~j ]C orr. > and estimate fc,j_ lil;corr .. from (|25|), 

(3) compare the values of <f>' q , k' jq _ 5 and <f> q . corr ., k jiq _j. corr , (j = 1, . . . , q - 1). If 
the relative difference is large, assume : 

'q; corr. > 
'3i1~3i corr. 

and repeat steps (2), (3) recursively until the relative difference between parameters 
before and after step (2) is small enough. Then go to the output, assuming 4> q := 4>' q , 
k — h' 

(OUTPUT) parameters </> 1; . . . , (f) q , k jt i (j + 1 = 1, . . . , q). 

Technical details and problems which appear when applying the algorithm to 
data were discussed in detail in . 

We have performed numerical simulations of the a—, p— and (p + a)— models || 
in order to test how the implementation algorithm works in practice. We generated 
10000 cascades of the 10 step length for the a— and (p+a)— models, and one cascade 
of the 10 step length for the p— model ^ for two different parameter sets separately. 

Implementation algorithm analized the data of the last cascade step. For each 
event it estimated the value of normalized intermittency exponents 4>2;norm., (p3;norm. 

), using the standard method (step 1) with recovering corrections 
included (steps 2,3). The selected results (for one set of parameters) are presented 
in Fig. 1 and in Tabs. 1, 2. 

For the a— model the histograms of 2; norm ., (p3- n orm. obtained in the standard 
method and the histograms with recovering corrections included are almost identical 
(see Figs, la, lb and Tab. 1). In this case the recovering corrections can be imple- 
mented better when one applies directly dedicated a— model recovering correction 
H i. e. if one substitutes coefficient kj A _j in (|TT|) by the product : kj q _j = 2^ ■2^ >q - j . 

3 It can be proved that for a given parameter set the p— model generates always the same values 
of the correlators and density moments 



k' 
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On the contrary, the implementation algorithm works well for the (p + a)— model 
(see Figs, lc, Id and Tab. 2). For the (p + a)— model the histogram with the 
recovering corrections included approximates well the theoretical value of normalized 
intermittency exponent. The histogram obtained by using the standard method is 
moved slightly to the left in comparison to the histogram with recovering corrections 
included. 

We have checked that for the p— model the theoretical values of normalized in- 
termittency exponents are estimated perfectly by both standard method and imple- 
mentation algorithm ||. 

It should be also mentioned that the histograms generated by the implementation 
algorithm (recovering corrections) are symmetric, in contrast to the standard ones, 
and their dispersions are of the same order as those derived for the standard method 
(cf. Tabs. 1, 2). 

To sum up we analyzed the estimation of intermittency exponents from the data 
which were generated by a multiplicative random cascading process. The following 
methods were applied: the standard method of cascade recovering (0) and the im- 
proved method which included recursively the recovering corrections. The improved 
method was applied in the form of the implementation algorithm. Numerical sim- 
ulations have been performed to check how both methods work in practice. The 
conclusions may be summarized as follows : 

(a) standard method of estimation of intermittency exponents does not apply 
for the whole class of multiplicative models : its accuracy depends on the specific 
properties of the model and its parameters. The method does not detect a conser- 
vation law if present in the model; 

(b) we propose an improved method of estimation of intermittency exponents. It 
exploits, in addition to the standard density (factorial) moments, the information on 
the bin-bin correlations, observed in the data and expressed in terms of the density 
correlators; 

(c) the method is formulated in the form of the recursive algorithm which, start- 
ing from the parameters obtained from the density moments and correlators, allows 
successive improvements of the result; 

(d) the method was tested in MC simulations which show that it is workable and 
indeed brings the experimental estimates closer to their theoretical values. Moreover, 
the improved distributions are symmetric with approximately the same dispersions 
as the uncorrected ones. 
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Figure captions 

Figs, la, lb Estimation of normalized intermittency exponents 4>2-,norm. and 03 ;norm . 
for a— model, using the standard method (dotted line), the improved method with 
the implementation algorithm (thin solid line), and dedicated a- corrections || 
(dashed line) compared with the theoretical values (solid line), performed for one 
set of a— model parameters : 

a\ = 0.8, a2 = 1.1, pi = 1/3 

Figs, lc, Id Estimation of normalized intermittency exponents <p2;norm. and (p3-norm. 
for (p + a)— model, using the standard method (dotted line), the improved method 
with the implementation algorithm (thin solid line), compared with the theoretical 
values (solid line), performed for one set of (p + a)— model parameters : 

a 2 i = 1 - a 2 i-i, P2i = V%-\ for % = 1, . . . , 10, 

ai = 0.2, as = 0.5, a 5 = 0.6, a? = 0.3, a 9 = 0.45, 

an = 0.25, ai 3 = 0.1, ai 5 = 0.15, an = 0.87, a w = 0.66, 

2pi = 0.05, 2p 3 = 0.15, 2p 5 = 0.25, 2p 7 = 0.40, 2p 9 = 0.05, 

2 Pll = 0.05, 2p 13 = 0.02, 2p 15 = 0.02, 2p 17 = 0.005, 2p 19 = 0.005, 
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Tab. 1 Estimation of normalized intermittency exponents 4>2-,norm. and 03 ;norm . and 
their dispersions for the a— model, using the standard method (second column), 
the improved method with the implementation algorithm (third column), and dedi- 
cated a— corrections || (fourth column), compared with the theoretical values (first 
column), performed for one set of a— model parameters (cf. Figs, la, lb) : 





theor. 


standard 


algorithm 


a-corr. 


02; norm. 


0.0285 


0.0251 ±0.004 


0.0246 ± 0.0033 


0.0288 ±0.004 


03; norm. 


0.0813 


0.0757 ±0.010 


0.0727 ±0.009 


0.0798 ±0.0111 



Tab. 2 Estimation of normalized intermittency exponents 4>2-norm. and 4>^ n orm. and 
their dispersions for the (p+a)— model, using the standard method (second column), 
the improved method with the implementation algorithm (third column), compared 
with the theoretical values (first column), performed for one set of (p ± a)— model 
parameters (cf. Figs, lc, Id) : 





theor. 


standard 


algorithm 


02: norm. 


0.177 


0.170 ±0.023 


0.173 ±0.029 


03; norm. 


0.478 


0.438 ±0.069 


0.470 ± 0.092 
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